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Voltage  instability  in  power  distribution  systems  can  result  in  voltage 
collapse  throughout  the  grid.  Today,  with  the  advanced  of  power  generation 
technology  from  renewable  sources,  concerns  of  utility  companies  are  much 
being  focused  on  the  stability  of  the  grid  when  there  is  an  integration  of 
distributed  generation  (DG)  in  the  system.  This  paper  presents  a  study  on  DG 
units’  placement  and  sizing  in  a  radial  distribution  network  by  using  a  pre¬ 
developed  index  called  Voltage  Stability  Condition  Index  (VSCI).  In  this 
paper,  VSCI  is  used  to  determine  DG  placement  candidates,  while  the  value 
of  power  losses  is  used  to  identify  the  best  DG  placement.  The  proposed 
method  is  tested  on  a  standard  33-bus  radial  distribution  network  and 
compared  with  existing  Ettehadi  and  Aman  methods.  The  effectiveness  of  the 
method  is  presented  in  terms  of  reduction  in  power  system  losses, 
maximization  of  system  loadability  and  voltage  quality  improvement.  Results 
show  that  VSCI  can  be  utilized  as  the  voltage  stability  indicator  for  DG 
placement  in  radial  distribution  power  system.  The  integration  of  DG  is 
found  to  improve  voltage  stability  by  increasing  the  system  loadability  and 
reducing  the  power  losses  of  the  network. 
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1.  INTRODUCTION 

Distributed  generation  attributes  to  generation  of  power  at  the  point  of  consumption  or  the  location 
of  load  [1,  2].  Distributed  generation  also  refers  to  small  scale  generation  that  is  not  directly  connected  to  the 
bulk  distribution  system  and  non-centrally  dispatched  [3].  In  published  works,  there  are  a  number  of  different 
terms  and  definitions  are  used  for  the  description  of  distributed  generation  [4].  The  term  ‘embedded 
generation’  is  often  used  in  Anglo-American  countries,  ‘dispersed  generation’  is  used  by  North  American 
countries  and  the  term  ‘decentralised  generation’  is  applied  in  Europe  and  parts  of  Asia  [4,  5].  The  name 
‘distributed  generation  (DG)’  is  recommended  to  be  used  all  over  the  world  based  on  the  literature  survey  [4, 
5].  Power  generation  at  the  location  of  load  decreases  the  losses  due  to  transmission  over  long  distances  as 
well  as  other  costs  associated  with  the  distribution  systems.  The  complexity  and  the  inefficiencies  from  the 
systems  are  also  minimized  with  the  emergence  of  DG  [6].  In  many  cases,  DG  can  supply  electricity  at  a 
lower  cost  with  improved  power  reliability  and  security  and  resulted  to  fewer  environmental  issues  compared 
to  conventional  generators  [7]. 

Voltage  stability  is  defined  as  the  ability  of  a  power  system  to  maintain  acceptable  voltage  levels  at 
all  buses  within  the  system  [8].  The  definition  includes  the  voltage  of  buses  under  normal  condition  and  after 
disturbance  in  the  network.  A  state  of  voltage  instability  is  experienced  in  a  system  when  there  is  a 
continuous  or  uncontrollable  drop  in  the  magnitude  of  voltage  either  after  a  disturbance,  load  increase  or  after 
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a  change  occurs  in  operating  conditions.  Voltage  instability  can  result  in  voltage  collapse  throughout  the  grid, 
in  which  case  the  voltage  falls  to  an  unacceptable,  low  value.  Voltage  collapse  is  a  common  occurrence  in 
power  systems  which  are  heavily  loaded,  faulted  or  having  shortage  in  reactive  power  [9,  10].  The  change  in 
voltage  is  directly  proportional  to  change  in  load  [11]. 

The  relation  of  grid  voltage  stability  and  grid  connected  DG  can  be  analysed  in  a  multitude  of 
methods.  As  much  as  the  reactive  power  support  sourced  from  capacitor  banks  can  be  utilised  to  reduce 
power  loss,  incorrect  placement  of  the  capacitors  can  further  increase  the  line  losses  [12].  Capacitors  are  also 
unable  to  support  the  network  with  real  power  when  the  situation  arises.  DG  systems  offer  a  better  option  as 
it  can  be  utilised  to  support  the  grid  with  both  real  and  reactive  power,  as  long  as  the  placement  and  the  sizing 
can  be  optimally  determined.  Inappropriate  sizing  and  location  of  distributed  generation  units  may  also 
increase  the  system  losses  as  well  as  the  operating  costs.  Therefore,  it  is  necessary  to  investigate  the  impacts 
of  DG  placement  and  sizing  to  the  system. 

Nannapaneni  et  al.  provide  a  comprehensive  study  to  evaluate  the  effectiveness  of  various  PQ  and 
QV  voltage  stability  indices  on  utilizing  distributed  generation  (DG)  units  optimally  within  certain  limits  and 
constraints  [13].  The  indices  are  used  to  identify  the  weakest  buses  and  calculate  the  stability  margin  of  a 
system.  The  study  concludes  that  QV  sensitivity  indices  need  to  be  considered  along  with  PV  sensitivity 
indices  for  the  allocation  of  DGs.  The  study  also  shows  how  the  selection  of  the  appropriate  voltage  stability 
index  varies  depending  on  the  DG  bus  model  in  the  simulation,  either  as  a  PV  or  PQ  model.  The  study 
proposes  a  more  direct  method  in  suggesting  QV  index  in  analysing  voltage  stability  as  QV  curve  method 
does  not  require  the  system  to  be  represented  as  a  two-bus  system.  However,  Q-V  curve  could  only  be  used 
for  a  system  with  a  constant  power  level,  thus  it  needs  P-V  curve  method  to  be  performed  as  well. 

In  a  research  by  Nahilia  and  Boudour,  a  static  power  flow  study  based  on  the  analysis  of  Djanet 
power  system  behaviour  in  terms  of  DG  contribution  and  voltage  profile  evolution  in  Algeria  [14].  A  new 
voltage  stability  index  is  used  to  identify  the  weakest  node  in  the  system.  This  paper  states  that  the  optimal 
DG  power  factor  must  be  closer  to  the  power  factor  of  the  combined  load  contribution  in  order  to  minimize 
power  losses.  The  results  of  this  research  show  that  the  integration  of  DG  improves  the  voltage  profile  and 
decreases  the  total  power  loss  in  the  system  significantly.  The  research,  however,  does  not  consider  the  sizing 
and  the  location  of  DG  in  the  study. 

Wang  et  al.  investigate  the  impact  of  wind  power  on  static  voltage  stability  [15].  The  probabilistic 
static  voltage  stability  of  power  systems  is  analysed  based  on  Probabilistic  Load  Flow  (PLF),  which  is  able  to 
take  uncertain  factors  into  consideration,  for  example  load  distribution  and  randomness  of  input  variables. 
Quasi-Monte  Carlo  (QMC)  simulation  using  Sobol  sequence  is  applied  to  solve  PLF  after  considering  the 
random  input  variables,  while  Nataf  transform  is  used  to  manage  the  correlation  of  wind  power.  The  results 
presented  in  this  paper  are  compared  with  the  data  generated  from  Monte-Carlo  Simulation  (MCS)  combined 
with  Simple  Random  Sampling  method  to  show  the  improvement  made.  However,  the  Nataf  transform  and  a 
Sobol  sequence  which  are  needed  for  this  method  makes  the  process  to  be  quite  tedious. 

Aman  et  al.  Has  introduced  a  new  index  for  the  placement  and  sizing  of  DG  units  [16].  The 
allocation  and  sizing  of  DG  units  in  this  research  using  the  proposed  index  is  to  determine  the  bus  with  the 
most  voltage  sensitivity  and  minimum  total  power  losses.  The  index,  power  stability  index  (PSI),  is 
developed  by  taking  stable  node  voltage  in  consideration.  The  simulation  results  show  that  the  DG  units 
improve  the  voltage  profile  of  the  network  and  increase  the  system  voltage  stability.  The  overall  system 
capacity  also  increases  while  reducing  the  line  losses.  The  proposed  algorithm  is  found  to  consume  less 
computation  time  compared  to  the  Golden  Section  Search  (GSS)  algorithm  by  50-60%.  However,  the 
performance  of  the  technique  depend  heavily  on  the  initial  guess  and  the  algorithm  is  developed  with  several 
theoretical  assumptions  such  as  continuity,  convexity  and  differentiability  which  may  not  be  suitable  for 
optimization  problems. 

Based  on  the  literature  reviews,  it  can  be  concluded  that  there  are  many  different  indices  used  for 
voltage  stability  analysis  as  well  as  many  different  approaches  for  DG  placement  and  sizing.  Some  of  the 
voltage  stability  indices  are  not  entirely  reliable  on  their  own  and  need  to  be  calculated  together  with  another 
index  for  it  to  be  accurate.  While  most  indices  are  comparatively  reliable,  some  of  the  approaches  pursued  for 
the  placement  and  sizing  of  DG  are  too  resource  extensive.  A  simpler  and  less  resource  extensive  method  is 
needed  for  the  voltage  stability  analysis  in  the  presence  of  DG  units.  This  paper  presents  the  ability  of 
Voltage  Stability  Condition  Indicator  (VSCI)  in  DG  placement.  Previously,  the  capability  of  VSCI  as  the 
voltage  stability  indicator  in  transmission  system  has  been  proven  in  [17-20].  In  this  paper,  VSCI  is  evaluated 
in  a  radial  distribution  system  network.  The  objective  of  DG  integration  in  this  paper  is  to  improve  the 
system  loadability  and  reduce  power  system  losses. 
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2.  RESEARCH  METHOD 

2.1.  Voltage  stability  condition  indicator  (VSCI) 

VSCI  has  been  proposed  in  [17]  and  is  derived  from  the  current  equation  in  a  complex  form  of  a 
simple  2-bus  Thevenin  equivalent  system.  For  a  typical  transmission  line,  VSCI  can  be  calculated  by 


VSCI  = 


AZn\p^+Q^) 


2x2 


<1 


(1) 


Where, 

P  and  Q  =  active  and  reactive  power  load, 

EjiiZ  S  =  Thevenin  voltage, 

ZjfiZ  p  =  Thevenin  impedance  or  the  total  equivalent  impedance  across  the  bus, 
and  d  and  P  are  their  corresponding  phase  angles. 

In  order  to  maintain  the  stability  of  the  system,  VSCI  value  must  be  less  than  or  equal  to  one  as  in 
equation  (1).  Any  bus  that  has  VSCI  value  near  to  one  implies  that  the  bus  is  approaching  its  collapse  point. 
Previously,  VSCI  has  been  tested  to  be  able  to  evaluate  the  voltage  stability  in  transmission  line  [17-20].  This 
paper  aims  to  test  the  capability  of  VSCI  to  evaluate  the  voltage  stability  in  radial  distribution  system. 


2.2.  DG  placement  and  sizing  based  on  VSI 

Today,  DGs  are  used  to  improve  the  traditional  electric  power  systems  [21-22].  These  technologies 
include  small  or  micro  gas  turbines,  fuel  cells,  wind  and  solar  energy,  biomass,  small  hydro-power  and 
others.  In  order  to  achieve  the  benefits  of  integrating  DG  units  into  the  power  system,  an  appropriate  place 
and  size  of  the  DG  units  installed  must  be  selected.  Installing  the  units  at  inappropriate  places  may  weaken 
the  stability  and  reliability  of  the  power  system  [23].  In  this  paper,  the  study  is  being  divided  into  two  phase. 
In  the  first  phase,  the  candidate  locations  of  DG  units  is  chosen  according  to  VSCI  value.  This  phase  is  very 
useful  to  minimize  the  size  to  be  used  in  the  optimization  process  for  the  second  phase  especially  for  very 
large  size  of  power  system  network.  Secondly,  the  best  location  and  sizing  of  the  DG  units  are  then 
determined  based  on  the  objective  function  which  are  in  this  paper  the  aim  is  to  minimize  the  power  losses 
and  maximize  the  system  loadability.  Figure  1  (a)  shows  the  flowchart  for  the  first  phase  of  the  study. 

In  the  second  phase  of  the  study  as  shown  in  Figure  1  (b),  the  ranking  of  the  weakest  buses  based  on 
the  VSCI  values  are  used  to  determine  the  best  placement  of  DG  units  to  be  installed  into  the  system.  The  top 
10  weak  buses  with  VSCI  above  0.95  are  used  to  place  a  DG  unit  at  each  location  and  the  effects  on  the 
network  power  losses  are  studied.  The  bus  that  give  the  minimum  value  of  power  losses  are  chosen  to  be  the 
optimal  location.  Next,  to  calculate  the  size  of  DG  unit  to  be  installed,  several  parameters  are  taken  into 
account  which  are  the  DG  penetration  percentage  (Pn),  power  factor  (pf)  value  of  the  system  and  number  of 
DG  units  (Num).  In  this  paper,  the  size  of  DG  are  determined  by  varying  the  penetration  percentage  and 
power  factor.  As  a  result,  in  this  paper,  the  best  penetration  level  of  DG  and  power  factor  value  of  the  system 
can  be  identified.  The  types  of  DG  used  are  made  constant  which  can  inject  both,  active  and  reactive  power. 
The  value  of  apparent  power,  S,  real  power,  P  and  reactive  power,  Q  to  be  supplied  by  the  DG  unit  is 
calculated  using  the  steps  below; 

a.  The  penetration  level,  Pn  and  the  total  apparent  load,  in  the  network  is  used  to  calculate  the 
apparent  power  to  be  supplied  by  the  DG  unit,  S^q  as  shown  in  the  formula  below; 

Pn  =  X  100%  (2) 

^load 


Pn  X 
100% 


(3) 


b.  The  real  power  to  be  supplied  by  the  DG  unit,  P^g  is  calculated  by  using  the  value  of  DG  apparent 
power  and  the  power  factor  set  for  the  DG  unit; 


Pdg  —  pf  ^  ^DG 


(4) 


c.  The  reactive  power  to  be  supplied  by  the  DG  unit,  is  calculated  using  the  formula  below; 


Qdg 


(5) 
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d.  The  real  and  reactive  power  to  be  supplied  by  the  DG  umt(s)  is  divided  equally  among  the  number  of 
units  to  be  integrated  in  the  network; 


p  = 

Num 


(6) 


Q  =  ^ 

Num 


(7) 


(a) 


(b) 


Figure  1 .  Flowchart  of  the  developed  algorithm  for;  (a)  DG  placement,  (b)  DG  sizing. 


3.  RESULTS  AND  DISCUSSION 

The  simulation  in  this  paper  was  computed  using  the  3  3 -bus  radial  distribution  network  [24].  This  is 
a  radial  distribution  system  with  33  buses,  5  tie  lines  or  looping  branches.  It  consists  of  1  generator  bus  as  a 
swing  bus,  32  load  buses  and  37  interconnecting  branches.  Table  1  shows  the  system  details  of  33-bus  radial 
distribution  network. 


Table  1.  33-bus  radial  distribution  system’s  details 


Total  Active  Load  in  kW 

Total  Reactive  Load  in 

Total  Apparent  Load  in 

System  Power 

(P) 

kVAR  (Q) 

kVA  (S) 

Factor  (pf) 

3715 

2300 

4369.35 

0.8502 

3.1.  Weak  buses  based  on  VSCI 

In  the  first  phase  of  research,  in  order  to  determine  weak  busses,  load  for  every  bussed  were 
increased  until  the  system  reached  its  maximum  loadability  limit.  VSCI  value  was  then  calculated  and  the 
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system  was  ranked  according  to  VSCI  value.  The  VSCI  values  that  close  to  unity  were  considered  as  weak 
buses.  Table  2  shows  the  rank  of  top  10  weak  busses  at  the  maximum  loadability  of  the  system.  The  results 
are  consistent  with  [25]  and  the  voltage  buses  magnitudes  which  give  almost  the  same  weak  buses  ranking. 


Table  2.  Top  10  busses  with  highest  VSCI  for  maximum  loadability 


Bus  No 

Voltage  (p.u.) 

VSCI 

18 

0.429 

0.9999 

17 

0.4337 

0.9849 

16 

0.4494 

0.9527 

32 

0.4609 

0.9443 

33 

0.4589 

0.9331 

31 

0.4676 

0.9265 

15 

0.4598 

0.9126 

14 

0.4705 

0.8828 

30 

0.4975 

0.8750 

13 

0.4878 

0.8646 

Real  Power  Losses  (MW) 

7.8640 

Reactive  Power  Losses  (MVAR) 

5.3385 

3.2.  DG  placement  and  sizing 

A  DG  unit  supplying  both  real  and  reactive  power  was  placed  at  each  weak  buses  ranked  as  the  top 
ten  based  on  VSCI.  In  this  paper,  the  type  of  DG  units  placed  in  the  network  can  inject  both  active  and 
reactive  power.  The  performance  of  DG  units  with  different  operating  power  factors  are  evaluated  based  on 
the  resulting  maximum  loadability  of  the  network,  which  are  shown  in  Table  3.  The  calculation  to  determine 
amount  of  DG  reactive  and  active  power  injected  is  as  discussed  in  Section  2.2.  In  the  beginning,  the  size  was 
set  to  40%  penetration  level  (Pn)  and  at  operating  power  factor  of  0.85  [25-26].  The  network  loadability  was 
set  constant  at  3.5,  which  was  the  maximum  loadability  of  the  base  case  network.  The  real  and  reactive 
power  losses  of  the  simulation  conducted  were  analysed.  From  Figure  3,  the  system  has  the  lowest  power 
losses  when  a  DG  unit  was  placed  at  buses  32,  33  and  31,  respectively.  Therefore,  in  this  paper,  buses  32,  33 
and  3 1  are  determined  as  the  best  location  for  the  DG  placement.  This  result  is  consistent  with  the  findings 
obtained  in  [25-26]. 

The  size  of  DG  was  determined  by  varying  the  penetration  level  (Pn)  of  DG.  In  this  paper,  the 
penetration  level  of  DG  was  varied  ranging  from  20%  to  80%,  from  one  unit  of  DG  up  to  three  units  of  DG 
evaluated  in  the  system.  Based  on  Table  3,  it  can  be  seen  that  the  network  has  the  highest  maximum 
loadability  when  the  DG  unit(s)  were  operating  at  0.85  power  factor.  This  is  due  to  the  DG  operating  power 
factor  is  being  closer  to  the  system  power  factor,  which  is  about  0.85  as  tabulated  in  Table  1.  This  result  is 
consistent  with  the  findings  reported  in  [26]  which  states  that  the  maximum  loadability  of  the  system  when 
being  integrated  with  DG  can  be  achieved  when  the  DG  operating  power  factor  is  similar  to  the  system 
power  factor.  Table  3  also  shows  that  although  the  size  and  the  number  of  DG  units  are  increased,  the 
maximum  loadability  achieved  are  more  less  the  same  for  placement  of  two  DG  units.  Hence,  considering  the 
installation  and  maintenance  cost,  two  DG  units  with  40%  penetration  level  are  recommended  for  33-bus 
radial  distribution  network.  This  result  is  again  consistent  with  the  finding  obtained  in  [26]  which  confirms 
the  capability  of  VSCI  and  the  value  of  power  losses  obtained  in  DG  placement  and  sizing. 


Table  3.  Comparison  of  maximum  loadability  of  different  DG  penetration  level  and  power  factor. 
DGs  =  1  DGs  =  2  DGs  =  3 

pf  DG  Penetration  (%)  DG  Penetration  (%)  DG  Penetration  (%) 


20 

40 

60 

80 

20 

40 

60 

80 

20 

40 

60 

80 

1 

3.89 

4.06 

4.19 

4.3 

3.95 

4.25 

4.52 

4.76 

3.94 

4.24 

4.52 

4.77 

0.95 

3.92 

4.12 

4.28 

4.42 

4.01 

4.38 

4.72 

5.04 

4 

4.37 

4.73 

5.07 

0.9 

3.93 

4.13 

4.3 

4.44 

4.02 

4.4 

4.77 

5.11 

4.02 

4.4 

4.77 

5.13 

0.85 

P.93 

4.14 

4.3 

4.44 

4.03 

4.42 

4.78 

5.14 

4.02 

4.41 

4.79 

5.16 

0.8 

3.93 

4.13 

4.3 

4.44 

4.03 

4.42 

4.78 

5.14 

4.02 

4.41 

4.79 

5.16 

0.75 

3.92 

4.13 

4.29 

4.43 

4.02 

4.41 

4.78 

5.13 

4.02 

4.41 

4.79 

5.15 
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Bus  Number 
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Figure  3.  Comparison  of  power  losses  for  a  unit  DG  placed  at  weak  buses. 


4.  CONCLUSION 

In  this  paper,  an  evaluation  of  a  voltage  stability  index  for  DG  placement  in  radial  distribution 
network  is  presented.  VSCI  is  used  to  determine  DG  placement  candidates,  while  power  losses  is  used  to 
choose  the  best  DG  places.  The  weak  buses  ranking  based  on  VSCI  gives  almost  similar  result  when  compare 
to  voltage  buses  magnitudes  and  also  findings  in  [25].  Therefore,  it  can  be  concluded  that  VSCI  can  be 
utilized  as  the  stability  indicator  of  radial  distribution  power  system.  The  study  is  conducted  on  the  33-bus 
radial  distribution  network,  and  the  results  show  the  impacts  of  DG,  both  in  total  power  loss  reduction  and 
maximum  loadability  improvement  under  stressed  operation.  Results  show  that,  with  40%  DG  penetration 
only,  power  losses  can  be  reduced  up  to  48%  and  maximum  loadability  of  the  system  can  increase  up  to  30%. 
From  above  results,  it  can  also  be  observed  that  the  performance  of  different  number  of  DG  units’  placement 
can  give  approximately  the  same  results.  This  result  is  again  consistent  with  the  findings  obtained  in  [26]. 
Therefore,  based  on  results  presented  in  this  paper,  VSCI  is  proven  to  be  able  to  determine  the  weakest  bus  in 
radial  distribution  system  for  DG  placement  and  sizing. 
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